Resolution of the B — > 7r7r, 7tX puzzles 
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We show that there exist uncanceled soft divergences in the kr factorization for nonfactorizable 
amplitudes of two-body nonleptonic B meson decays, similar to those identified in hadron hadropro- 
duction. These divergences can be grouped into a soft factor using the eikonal approximation, which 
is then treated as an additional nonperturbative input in the perturbative QCD formalism. View- 
ing the special role of the pion as a qq bound state and as a pseudo Nambu-Goldstone boson, we 
postulate that the soft effect associated with it is significant. This soft factor enhances the non- 
factorizable color-suppressed tree amplitudes, such that the branching ratios _B(7r°7r°) and B(-K°p°) 
are increased under the constraint of the B(p°p°) data, the difference between the direct CP asym- 
metries Acp(tv^ K ) and Acp(tt°K ) is enlarged, and the mixing-induced CP asymmetry S 7r a Ks 
is reduced. Namely, the known nn and -kK puzzles can be resolved simultaneously. 

PACS numbers: 13.25.Hw, 12.38.Bx, 12.39.St 



The more precise data of the B — > 7T7t, ttK decays 
have sharpened the discrepancies with the theoretical 
predictions from the factorization approaches, such as 
the perturbative QCD (PQCD) approach based on 
the kx factorization theorem [l|, 0]. The observed 
B° — > 7T 7r° branching ratio Q remains several times 
larger than the naive expectation. The direct CP 
asymmetry of the B^ — > decays differs dra- 

matically from that of B° — ► tt^K^ decays. There 
is a deviation between the extractions of the stan- 
dard model parameter sin(20i) from the penguin- 
dominated B° — > tt°Ks modes and from the tree- 
dominated b — > ccs modes. All these discrepancies are 
closely related to the color-suppressed tree amplitudes 
C 0. The B° -> vrV branching ratios from PQCD 
and QCD factorization (QCDF), being sensitive to C, 
are also much lower than the data p, Q . However, 
the estimate of C from PQCD is well consistent with 
the measured B° — > p°p° branching ratio 0] • Propos- 
als resorting to new physics Q mainly resolve the ttK 
puzzle without addressing the peculiar feature of C in 
the 7r°7r°, tt°p°, and p°p° modes, while those to QCD 
effects are usually strongly constrained by the pp data 
. It indicates the difficulty of resolving the 7T7t and 
%K puzzles simultaneously. 

The color-suppressed tree amplitude C seems to be 
an important but the least understood quantity in B 
meson decays. Viewing that all the puzzles appear 
in the C-sensitive quantities, we shall carefully inves- 
tigate QCD effects on C, and their impact on the 
B — > 7T7T, ttK decays. Once a mechanism identified 
for C respects the conventional factorization theorem, 
it is unlikely to be a resolution due to the B — > pp 
constraint mentioned above 0]. That is the reason 
the higher-order corrections calculated in QCDF , 
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FIG. 1: LO diagrams for a nonfactorizable amplitude. 



which obey the collinear factorization, cannot resolve 
the 7T7T puzzle. It has been pointed out by Collins 
and Qiu [TTf| that the kr factorization breaks down 
in complicated QCD processes like high-pT hadron 
hadroproduction. To factorize the collinear gluons as- 
sociated with, say, one of the initial-state hadrons, one 
needs to eikonalize the valence quark lines to which 
the collinear gluons attach. Those eikonal lines, i.e., 
Wilson lines from another initial-state hadron and the 
final-state hadrons, should cancel in order to have 
the universality of the considered parton distribution 
function. However, the required cancellation is not ex- 
act in the kx factorization, though it is in the collinear 
factorization. The kx factorization still holds for sim- 
ple processes like deeply inelastic scattering (DIS), 
which does not involve the Wilson lines from the other 
hadrons. 

The above observation provides a clue for resolv- 
ing the 7T7r and ttK puzzles. It is easy to see that 
a factorizable amplitude, involving only a B me- 
son transition form factor, mimics simple DIS, and 



a nonfactorizable 27 1 amplitude, involving dynamics 
of three hadrons, mimics the complicated hadron 
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FIG. 2: NLO diagrams for Fig. QJa) that are relevant to 
the factorization of the A'h meson wave function. 



hadroproduction. The kr factorization for a factoriz- 
able B meson decay amplitude has been proved [Hj]. 
The kx factorization for a nonfactorizable amplitude 
has not, though it has been widely employed in the 
PQCD analysis. Below we shall identify the resid- 
ual infrared divergence in the fcr factorization for 
a nonfactorizable amplitude at one loop. Contrary 
to high-pT hadron hadroproduction, this residual in- 
frared divergence can be factorized into a soft factor 
in two-body nonleptonic B meson decays, following 
the procedure in [l3| , so that the universality of a /in- 
dependent meson wave function is restored. A non- 
factorizable amplitude then remains calculable in the 
PQCD approach after parameterizing the soft factor. 

Consider the B{P B ) -> M 1 {P 1 )M 2 (P 2 ) decay, 
where Pb, 1,2 represent the momenta of the B, Mi, and 
M 2 mesons, respectively. Start with the leading-order 
(LO) nonfactorizable emission diagram in Fig. [TJa) 
resulting from the operator 2 We add a ra- 

diative gluon of momentum I emitted by the valence 
quark in M%. The attachment of the radiative gluon 
to the b quark line shown in Fig. Ufa) leads to a Wilson 
line from infinity to the origin, i.e, the weak vertex. 
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FIG. 3: Two-particle reduced NLO diagrams for the M2 
meson wave function. 



This piece is factorized in color flow by itself with the 
color factor Cf ■ The attachment to the hard gluon in 
Fig- Hlb) generates two Wilson lines, one of which runs 
from the position z 2 of the valence anti-quark in M 2 to 
infinity. The attachments to the virtual anti-quark in 
Fig- EJc) and to the valence quark in the Mi meson in 
Fig. [2jd) also generate the Wilson line running from 
z 2 to infinity. The combination of these three pieces 
with the same Wilson line is factorized in color flow. 
As to the next-to-leading-order (NLO) two-particle re- 
ducible diagrams, such as the self-energy correction to 
the valence quark in Fig.^a) and the gluon exchange 
between the valence quark and the valence anti-quark 
in Fig. [3jb), their factorization into the M 2 meson 
wave function is straightforward. 

The other attachments, such as Figs. [He) andO(f), 
should cancel in order to have the universality of the 
M 2 meson wave function. Eikonalizing the spectator 
quark line on the B meson side, Fig. [^e) gives the 
integral 
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where the vector n+ = (1, 0, Or) is along the momen- 
tum Pi, the ... denotes the rest of the integrand, and 
75^*2 comes from the twist-2 structure of the M 2 me- 
son wave function. Another piece from Fig. [2Kb) con- 
tains the Wilson line running from the position of the 
spectator in the Mi meson to infinity. This piece, 
together with Fig.f^f), leads to the loop integral the 
same as Eq. (JIJ , but with the eikonal propagator being 
replaced by 1/ (— n+ • l + ie). Employing the principal- 
value prescription 
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we find a NLO residual soft divergence from the 
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Glauber region with l~ — 0, which seems to violate 
the universality of the M 2 meson wave function. 

The residual soft divergence may not cause trou- 
ble, if the LO amplitude A^ ** is real. Expanding the 
decay width up to NLO, we have 



\M\ 2 = |7W (0) | 2 + 2Re[M iQ) M {1) *] 



(3) 



According to Eq. ([5]), the residual soft divergence will 
be purely imaginary, if Ai^ is real, so it does not 
survive in the second term Re[7W^ ^ M^*]. This is 
the reason the diagrams Fig. [5] are calculable in the 
collincar factorization, and result is mainly real in 
QCDF [ig. is complex in the factoriza- 

tion, because the partons carry transverse momenta, 
and the internal lines can go on mass shell at finite 
momentum fractions [Hj]. Then the residual soft di- 
vergence contributes to Re[M w M {1) *}. It has been 
known that the end-point singularity from a vanishing 
momentum fraction is developed at twist 3 in QCDF 
for a nonfactorizable decay amplitude. The presence 
of the end-point singularity signals the inappropriate 
theoretical framework, while the residual soft diver- 
gence in the kx factorization for a nonfactorizalbe B 
meson decay amplitude can be factorized and be han- 
dled by introducing an additional nonperturbative in- 
put as shown below. 

We eikonalize the valence quark line of the M 2 me- 
son according to the observation in [l3j : for P 2 in the 
minus direction, the poles 
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from the denominators (/c 2 + I) 2 + ie and (P 2 — &2 + 
k\ — I) 2 + ie in Eq. (JTJ) , respectively, are apart from 
each other by O(fc^). As long as k^ is greater than 
O(Aqcd), Aqcd being the QCD scale, we can deform 
the contour of l + , such that / + remains larger than 
O(Aqcd), and the approximation 

fc 2 -/+ > 0(m B A QCD ) » |k 2T + 1t| 2 ~ 0(A2j CD ), (6) 

holds. We then have the eikonalization of the valence 
quark line even in the Glauber region, 
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Picking up the pole l + = — ie from the above eikonal 
propagator, the residual soft divergence in Fig. [2] is 
factorized into the convolution 
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FIG. 4: A NLO diagram that does not contribute the resid- 
ual soft divergence. 



where Ai^ denotes the LO amplitude from Fig.QJa). 

If the radiative gluon is emitted by the valence anti- 
quark of the M 2 meson in Fig.[lja), the corresponding 
loop integral does not generate the residual soft diver- 
gence. Take the NLO diagram in Fig.2]as an example, 
where the radiative gluon attaches to the spectator of 
the Mi meson. All the poles of l + are located in the 
lower half plane, and not pinched, so the integral di- 
minishes. That is, these diagrams contribute only the 
infrared divergence that ought to be absorbed into the 
M 2 meson wave function. 

Adding the second radiative gluon emitted by the 
valence quark in the M 2 meson, we have 
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whose derivation is similar to the two-loop analysis 
in [llj . Extending the analysis to higher orders, and 

neglecting the dependence of on the transverse 

loop momenta It, the divergence can be summed into 
an exponential factor, I a ~ e tSc Ai^ with S e being a 
parameter. Considering the soft corrections to the LO 
diagram in Fig. \V[b) from the gluons radiated by the 



valence anti-quark in M 2 , we obtain If, 
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where the minus sign is attributed to the radiation 
from the anti-quark. The above modified fcy factor- 
ization formalism also applies to nonfactorizable am- 
plitudes from other tree and penguin operators, and 
to nonfactorizable annihilation amplitudes. 

The soft factor has a dynamical origin similar to 
that of a meson wave function: the former (latter) ab- 
sorbs the soft (collinear) gluons. The behavior of the 
soft factor can be obtained either by nonperturbative 
methods or from experimental data. We then study 
its impact on nonleptonic two-body B meson decays. 
The color-suppressed tree amplitudes are small at LO 
due to the small Wilson coefficient a 2 for the factor- 
izable contribution and to the pair cancellation be- 
tween Figs.[lla) and[ljb) for the nonfactorizable con- 
tribution. The presence of the soft factor can convert 
the destructive interference in Fig. [T] into a construc- 
tive one, leading to strong enhancement. Its effect on 
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the color- allowed tree and penguin (including annihi- 
lation) amplitudes should be minor, which are domi- 
nated by factorizable contributions. 

We seek experimental constraints on the soft factor 
by comparing the data [|[ (in units of 10~ 6 ) 



S(7r°7r°) 
B(pV) 



1.55 ±0.19, [(0.29 
2.0 ±0.5, [«o.r 
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0.74 
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[(0.921^°)], (10) 



with the NLO PQCD predictions in the square brack- 
ets, which are quoted from (ill. Jll\\. and 0, respec- 
tively. The results from QCDF [1,111] are similar. The 
above comparison motivates us to postulate that the 
soft effect is significant (negligible) in the decays with 
M2 being a pseudoscalar (vector) meson. That is, we 
associate a soft factor with M2 = tt, K, but not with 
Mi = p (the soft effect associated with the kaon is 
not crucial actually). Because the B° —* tt°p° decay 
involves both types of amplitudes with the pion and 
the p meson as M2, it is natural that the discrepancy 
is in between as indicated by Eq. (fTO)) . The larger soft 
effect from the multi-parton states in the pion than 
in the p meson can be understood by means of the 
simultaneous role of the pion as a qq bound state and 
as a Nambu-Goldstone (NG) boson the valence 
quark and anti-quark of the pion are separated by a 
short distance, like those of the p meson, in order to 
reduce the confinement potential energy. The multi- 
parton states of the pion spread over a huge space-time 
in order to meet the role of a massless NG boson. 

The factorization formulas for the B — > tttt, irp, ttK, 
and pK decays can be found in [l|| [2(| • According to 
our derivation, we multiply the b quark nonfactoriz- 
able emission amplitudes, both tree and penguin, by 
e lSe (e~ iSe ) with the hard gluon being emitted by the 
valence anti-quark (quark) in M.2- The dependence 
on S e of those C-sensitive quantities is displayed in 
Fig. |3 The branching ratios B(7r°7r°) and B(n°p°) 
grow quickly with decreasing S e from the NLO PQCD 
values in Eq. (fTOj) . and become close to the data 
when S e reaches — tt/2. The direct CP asymmetry 
A C p(7r X ± ) increases from the NLO PQCD result 
around —0.01 [l6[ to above 0.05 for S e < — rr/4, whose 
agreement with the data Acpi^K^) = 0.050 ±0.025 
[3| is satisfactory. The deviation of the mixing induced 



CP asymmetry AS n o Ks = S v o Ks 
the NLO PQCD value +0.07 to - 
Compared to the data S 7: o Ks 



- S C cs descends from 
0.04 for S e = -tt/2. 
0.57 ± 0.17 and 



Sees = 0.672 ± 0.024 the consistency has been im- 
proved. The ratio C/T = 0.53e~ 22i for the B -> tttt 
decays, corresponding to S e ~ —tt/2, is consistent 
with the extraction in [4|, T being the color-allowed 
tree amplitude. An equivalent viewpoint is that the 



B(tt°tt°) data constrain S e ~ —tt/2, which then leads 
to the predictions for other quantities in Fig. [5] 

We have confirmed that B(n Zf TT ± ), B(7r°7r ± ) and 
all B(ttK) change slightly from those in (l6| . since 
they are less sensitive to C. Acp{tt t K^) remains 
around —0.1 1] for arbitrary S e , and in agreement 
with the data Acp^K^ = -0.098±g;g|f [|. The 
small variation of the curve is attributed to the soft ef- 
fect on the nonfactorizable color-allowed and penguin 
contributions. S p a Ks does not change much, because 
M2 = p in this case, and the involved C is not mod- 
ified. The NLO PQCD prediction AS p o Ks « -0.15 
20] is consistent with the data S p o Ks = 0.631q 2i [jj. 
For those penguin-dominated two-body modes with- 
out involving C, like B — > 4>K, their mixing-induced 
CP asymmetries are not affected cither. 

Our resolution differs from those based on new 
physics models, such as the fourth-generation model 
|2ll . |22| . in which it is the electroweak penguin am- 
plitude that is enhanced. These proposals, in which 
new weak phases are introduced, change S p o^ s and 
S^Ks- O ur proposal, which does not involve multi- 
ple intermediate states, differs from the elastic rescat- 
tering models for final-state interaction [23]. A large 
C has been generated through the charge exchange 
mechanism in [24]]. Since the QCDF approach was em- 
ployed, the parameter scenario "S4" [6j or the inelas- 
tic scattering [25| has to be incorporated in order to 
get the correct result for Ac p(tt^ K^) . The exchange 
mechanism, which requires turning of two energetic 
quarks into opposite directions, is suppressed accord- 
ing to the factorization theorem. Moreover, AS n o Ks 
remains positive in [24| . 

In this letter we have identified the uncanceled soft 
divergences in the kx factorization for the nonfactor- 
izable B meson decay amplitudes, which are similar 
to those existent in hadron hadroproduction. The 
divergences are factorizable and demand the intro- 
duction of the additional soft factor, under which 
we have computed C in the PQCD approach, and 
found a simultaneous resolution of all the puzzles: 
B(tt°tt°) and B(tt°p°) are enhanced, the difference be- 
tween Acp{tt t K ± ) and Acp(tt°K ± ) is enlarged, and 
AS n a Ks is reduced for the single soft parameter S e 
around —tt/2. The constraint on C from the B — > pp 
data is evaded, because of the special role of the pion 
as a qq bound state and as a pseudo NG boson. The 
mechanism identified here can be verified or falsified 
by more precise data in the future. Based on our ob- 
servation, we agree that the B — > ttK data have not 
yet revealed a new physics signal [26j j . 
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Theoretical Sciences and National Science Council of 
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FIG. 5: S e dependence of (a) B(7r°7r°) and B(ir p ) in units of 10" 6 , of (b) Acp^ ^) and Acp^K*), and of (c) 
AS 7r o Ks and AS p o Kg . The data (horizontal bands) for AS are not shown due to their large errors. 
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